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Abstract. We review recent results on a mean-field model for rela- 
tivistic electrons in atoms and molecules, which allows to describe at 
the same time the self-consistent behavior of the polarized Dirac sea. 
We quickly derive this model from Quantum Electrodynamics and state 
the existence of solutions, imposing an ultraviolet cut-off A. We then 
discuss the limit A — 5* oo in detail, by resorting to charge renormaliza- 
tion. 
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For heavy atoms, it is necessary to take relativistic effects into account. 
However there is no equivalent of the well-known TV-body (non-relativistic) 
Schrodinger theory involving the Dirac operator, because of its negative 
spectrum. The correct theory is Quantum Electrodynamics (QED). This 
theory has a remarkable predictive power but its description in terms of 
perturbation theory restricts its range of applicability. In fact a mathemat- 
ically consistent formulation of the nonperturbative theory is still unknown. 
On the other hand, effective models deduced from nonrelativistic theories 
(like the Dirac-Hartree-Fock model |40|. 116]) suffer from inconsistencies: for 
instance a ground state never minimizes the physical energy which is always 
unbounded from below. 

Here we present an effective model based on a physical energy which can 
be minimized to obtain the ground state in a chosen charge sector. Our 
model describes the behavior of a finite number of particles (electrons), cou- 
pled to that of the Dirac sea which can become polarized. Our existence 
results are fully non-perturbative. Like in QED, the model contains diver- 
gences which have to be removed by renormalization. 

We review several results obtained in collaboration with Christian Hainzl, 
Philippe Gravejat, Eric Sere and Jan Philip Solovej. These works have 
already been summarized in [23] and in the fourth chapter of [29J, to which 
the interested reader is refered for more details. 

1. A NONLINEAR DlRAC EQUATION 

We present a mean-field model describing the self-consistent behavior of 
a finite number of 'real' electrons in an atom or a molecule, and, simul- 
taneously, of the infinitely many 'virtual' electrons of the Dirac sea. The 
state of the system is described by a one-body density matrix P which is 
a self-adjoint operator acting on the Hilbert space fj := L 2 (1R 2 ,C 4 ), and 
satisfying the constraint < P < 1 (1 denotes the identity operator on Sj). 
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The operator P describes the whole system consisting of the real and virtual 
electrons. We are interested in the following stationary equation [5} [2~5]: 

(1) 

In this section we explain the meaning of this equation at a formal level, 
before turning to rigorous results. 

The operator D is a mean-field Hamiltonian which is seen by all the 
particles. The first term 

D° :=a-(-iV) + P 

is the usual free Dirac operator |41] (with a = {a\, «2, a^) where ctj are the 
usual Dirac matrices). For the sake of simplicity we have chosen units in 
which the speed of light is c = 1 and the mass of the electron is m = 1. 
The operator D° satisfies (D ) 2 = —A + 1 and its spectrum is o~(D°) = 
(— oo, —1] U [l,oo). The second term in the formula of D is the Coulomb 
potential induced by both a fixed external density of charge v (typically 
v = Z5q for a pointwise nucleus of charge Z located at the origin), and the 
self-consistent density pp_i/2 of the electrons (defined below). The number 
a = e 2 which is the square of the (bare) charge e of the electron, is called 
the bare coupling constant. It will be renormalized later. 

The third term Xp in the definition of D is an exchange operator whose 
form depends on the chosen model. In Hartree-Fock (HF) theory [23] , we 
have 

,r , \ (P-l/2)(x,y) 

X P (x, y) = -a K - , 1 >\ ,U) (2) 

f - v\ 

which is called the exchange term. In Relativistic Density Functional Theory 

USE], 

Qp 

X P = -^-{pP-1/2) (3) 

is the derivative of a chosen effective exchange-correlation functional, which 
depends only on the density pp_i/2- In reduced Hartree-Fock (rHF) theory, 
we simply take 



X P = 0. 



For the sake of clarity, we will mainly present the mathematical results 
that have been obtained in the simplest case of Xp = and we will only 
make comments on the Hartree-Fock case ([2]). The exchange-correlation 
approximation ([3]) has not been considered rigorously so far. 
We use the notation X(-oc,/t) 

(D) to denote the spectral projector of D 
associated with the interval (— 00,/i). Hence Equation ([1]) means that the 
electrons of the system fill all the energies of the mean-field Hamiltonian 
D, up to the Fermi level \x S (—1,1). In practice we choose the chemical 
potential /i to fix the total charge of the systemQ We have added in (fl|) the 
possibility of having a density matrix < 5 < X{y,}(D) a ^ t ne Fermi level, as 



4f the external field is not too strong, fixing the charge is the same as fixing the number 
of electrons. However in strong fields, electron-positron pairs can be created and fixing 
the charge might not lead to the expected number of electrons. 
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is usually done in reduced Hartree-Fock theory [39] . So the operator P is not 
necessarily a projector but we still use the letter P for convenience. Later 
we will restrict ourselves to the case of P being an orthogonal projector. 

Equation ([1]) is well-known in the physical literature. A model of the 
same form (with an exchange term Xp different from ([2])) was proposed by 
Chaix and Iracane in [5]. Relativistic Density Functional Theory aims at 
solving the same Equation ([1]) with Xp given by ([3]) and additional classical 
electromagnetic terms accounting for the interactions with photons, see, e.g., 
[Po| Eq. (6.2)] and [HI Eq. (62)]. Dirac already considered in [10] the first 
order term obtained from (pQ) in an expansion in powers of a, assuming 
X P = 0. 

Let us now elaborate on the exact meaning of pp-\/2- The charge density 
of an operator A : S) — > S) with integral kernel A(x, y) a y is formally defined 
as pa{x) = Yla=l x)o,a = Tr C 4(A(x, x)). In usual Hartree-Fock theory, 
the charge density is pp(x). However, as there are infinitely many particles, 
this does not make sense here. In ([1]), the subtraction of half the identity 
is a convenient way to give a meaning to the density, independently of any 
reference, as we will explain later. One has formally, when P is a projector, 

Pp-i/2(x) = P p- P ± (x) = \Y^ l^ r ( x )! 2 - Wti*)? 

where {ip7}i>l is an orthonormal basis of PS) and {<f>t}i>l is an orthonormal 
basis of (1 — P)Sj. As was explained in [23] (see also Section [2]), subtract- 
ing 1/2 to the density matrix P renders the model invariant under charge 
conjugation. 

The free vacuum. When there is no external field (y = 0) and when 
Xp = 0, Equation ([1]) has an obvious solution for any p E (—1,1), the state 
made of all electrons with negative energjQ 

P = P° :=X(-oa,o)(D°), 
in accordance with Dirac's ideas [3 El [9]. Indeed p P o_ 1 / 2 = 0, as is seen by 
writing in the Fourier representation 

and using that the Dirac matrices are trace-less. This shows the usefulness 
of the subtraction of half the identity to P, since the free vacuum now 
has a vanishing density. 

For a general state P, we can use this to write (formally): 

Pp-i/2 = Pp-i/2 ~ PpO-1/2 = Pp-po- (4) 

When P belongs to a suitable class of perturbations of P° (for instance 
when P — P®_ is locally trace-class), the density p P _ P o is a well-defined 
mathematical object. We will give below natural conditions which garantee 
that P — P_ has a well-defined density in our context. 



2 In the Hartree-Fock case @, the free Dirac sea P = P° is no more a solution of 
|T]) when v = 0. The Hartree-Fock free vacuum solving the nonlinear equation (jTJ) was 
constructed in [311124] , assuming an ultraviolet cut-off. 
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Figure 1. State of the system depending on the value of the 
chemical potential fx. 



Electrons interacting with the polarized vacuum. With external field 
(v ^ 0), Equation ([1]) models a system of electrons in the presence of a 
nucleus and with a self-consistent polarized Dirac sea. The number of 'real' 
electrons in the system will depend on the value of fi. Typically (for not 
too strong fields) when \x = 0, one obtains the ground state of the polarized 
vacuum in the presence of u, without any real electron. On the other hand 
fx > in general leads to systems with a finite number of real electrons (Fig. 

ID. 

Note that for a generic operator < P < 1 there is no natural distinction 
between real and virtual electrons. It is only for a solution of Equation 
(JT]) that we can interpret the eigenfunctions corresponding to the positive 
eigenvalues of D as describing 'real' electrons, and the rest of the spectrum 
as being the Dirac sea. 

When fi > (and 5 = 0), the ./V filled eigenfunctions of D corresponding 
to the eigenvalues in [0, fi) solve the following system of nonlinear equations: 



x ' x 



for i = 1,...,N. This equation has the same form as the well-known Dirac- 
Hartree-Fock equations [lQl[16], without exchange term, and with an addi- 
tional vacuum polarization potential induced by the vacuum 
X(-oo,o)0^)- This remark was used by Chaix and Iracane in [5] as a justifi- 
cation to the Dirac-Hartree-Fock model. 

Most of the material of this section is purely formal and many objects (like 
PP—1/2) do not really make sense as such. In Section [3] we will introduce an 
ultraviolet cut-off A and present rigorous results. We however first explain 
how the formal Equation (fTJ) can be derived from Quantum Electrodynamics 
(QED). 

2. Derivation from Quantum Electrodynamics 

In this section we derive Equation ([1]) from first principles. We start with 
the formal QED Hamiltonian written in Coulomb gauge, in the presence of 
an external density of charge v and an external magnetic potential a, see 
[251 EH M EH SJ: 

«- = /*»(«■ {-*V-^(A W + „ W}+m /0*( 



a 



x x) dx 

p P^dx d y + U[ P 4^ d xdy + H f . (5) 
\x-y\ 2JJ \x-y\ 
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Here ^(x) is the second quantized field operator which annihilates an elec- 
tron at x and satisfies the anticommutation relation 

V*{x) V(y) v + *(y) v V*{x) a = 2S aiV S(x - y). (6) 

In the formula of M u ' a , p(x) is the density operator defined by 

P(x) = E ln(x) ;*' wl m 

where [a, b] = ab — ba. The operator Hf describes the kinetic energy of the 
photons: 

H f = — [ (|V x A(x)\ 2 + \E t (x)\ 2 ) dx = V / dk\k\a\{k)a x {k) + Cst 

(Cst indicates a constant which diverges in infinite volume). The operators 
A{x) and Et{x) are the electromagnetic field operators for the photons and 
a* x (k) is the creation operator of a photon with momentum k and polarization 
A. The Hamiltonian M u ' a formally acts on the Fock space J- = J- c ® J-" p h 
where J- e is the fermionic Fock space for the electrons and J-ph is the bosonic 
Fock space for the photons. 

We emphasize that ([5]) does not contain any normal-ordering or notion 
of (bare) electrons and positrons: ^f(x) can annihilate electrons of negative 
kinetic energy. The distinction between electrons and positrons should be a 
result of the theory and not an input. The commutator used in the formula 
([7]) of p(x) is a kind of renormalization, independent of any reference. It 
is due to Heisenberg [25] (see also [331 Eq. (96)]) and it is necessary for 
a covariant formulation of QED, see [33, Eq. (1.14)] and [IH Eq. (38)]. 
More precisely, the Hamiltonian M u ' a possesses the interesting property of 
being invariant under charge conjugation since the following relations hold 
formally 

Cp(x)C" 1 = -p(x), CW'X- 1 = M~ u ' a , 
where C is the charge conjugation operator acting on the Fock space. 

We now make two approximations: (i) we neglect photons and assume 
there is no external magnetic field, a = 0; (ii) we work in a mean-field 
theory, i.e. we restrict the Hamiltonian M u,a to (generalized) Hartree-Fock 
states. 

Let us recall that the electronic one-body density matrix (two point func- 
tion) of any electronic state |0) G 7 e is defined as 

and it satisfies < P < 1. Generalized Hartree-Fock states form a subset of 
(mixed) states which are completely determined by their density matrix P, 
see [3]. The value of any product of creation and annihilation operators is 
computed by means of Wick's formula. The energy of a Hartree-Fock state 
|HF) (g> |0) (with |0) G .F p h being the photonic vacuum) is [24"] 

=£gp(P-l/2) + Cst 
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where Cst is a constant (diverging in the infinite volume limit) and 
££ p (P-l/2) =Tr(I> (P- l/2))-a If Pp -y 2 ^^ dxdy 



2JJ \x-y\ 2JJ \x-y\ 



a 



The reader can recognize in (J8J) the well-known Hartree-Fock energy 
3j, but applied to the "renormalized" density matrix P — 1/2 instead of 
the usual density matrix P. The last two terms of the first line are respec- 
tively the kinetic energy and the interaction energy of the electrons with 
the external potential induced by the charge distribution v. In the second 
line appear respectively the so-called direct and exchange terms. In Rela- 
tivistic Density Functional Theory (RDFT) |144 1 15 j . the exchange term is 
approximated by an exchange-correlation functional F^ c {pp -1/2) whereas in 
reduced Hartree-Fock theory, the exchange term is simply dropped. 

Writing the first and second order stationarity conditions with respect to 
the density matrix P leads to the nonlinear equation ([1]) with [i = 0. The 
equation with fj, ^ is obtained by replacing D° by D° — fj,. Again our 
derivation is formal but (in the Hartree-Fock case) this was made rigorous 
by means of a thermodynamic limit in 



Remark 1. Instead of the vacuum, one can take a coherent state for the 
photons. This leads to a classical unknown magnetic field A{x) interacting 
with the particles. So far, there are no mathematical results on such a model. 



3. Existence and non existence of solutions 

In the presence of an external field (y 7^ 0), Equation ([T]) has no solu- 
tion in any 'reasonable' Banach space [21] and it is necessary to introduce 
an ultraviolet regularization parameter A. The simplest method (although 
probably not optimal regarding regularity issues [IT]) is to impose a cut-off 
at the level of the Hilbert space, that is to replace fj = L 2 (M 3 ; C ) by 

:= {/ G L 2 (M 3 ;C 4 ), supp(/) C B(0;A)} 
and to solve, instead of (HJ, the regularized equation in ft a'- 

P = X(-oo,m) (D) + S 

D = U A (l>° + a(p P _ P o - v) * Ix]" 1 ) n A 

where Ha is the orthogonal projector onto f)A in £)• We take Xp = in the 
rest of the paper, that is we work in the reduced Hartree-Fock approxima- 
tion. Note that we have used dU to replace pp-x/2 by p P _ P o. 

Existence of solutions. Existence of solutions to ([9]) was proved in |21] 
for /i = and in [T7] for fi G (—1,1), for all values of the coupling constant 
a > 0. The precise statement of this nonperturbative result is the following: 

Theorem 1 (Nonperturbative existence of solutions to ©, [HIE]). As- 
sume that a > 0, A > and pi £ (—1,1) are given. Let v he in the so-called 



(9) 



RENORMALIZATION OF DIRAC'S POLARIZED VACUUM 



7 



Coulomb space: 

C:=if : ^J£T 2 |/(*0| 2 ^<°o}. 

Then, Equation ([9]) has at least one solution P such that 
P-F^6 2 (%), P£(P-P°)P£ G 6i(£ A ), Pp^po G CfU 2 (IR 3 ). (10) 
All such solutions share the same density Pp_po ■ 

In (fTU|) . and @2(-?)a) are respectively the spaces of trace-class and 

Hilbert-Schmidt operators [38] on f)A, and P? = 1 — P°. Note that thanks 
to the uniqueness of Pp_po , the mean-field operator Z) is also unique and 

only 5 can differ between two solutions of ([9]). 

Let us mention that it is natural to look for a solution of ([9]) such that 
P — P° is a Hilbert-Schmidt operator on 53a- If P is a projector, the Shale- 
Stinespring theorem [37] tells us that P yields a Fock representation equiv- 
alent to that of P2_. Even when P is not a projector, it will be associated 
with a unique Bogoliubov mixed state in the Fock space representation of 

. This is a mathematical formulation of the statement that P should not 
be too far from P°. 

The method used in |2H ITT] to prove Theorem [T] was to identify solutions 
of © with minimizers of the so-called reduced Bogoliubov-Dirac-Fock energy 
which is nothing but the formal difference between the reduced Hartree-Fock 
energy of P and that of the reference state P°. A formal calculation yields 
with Q = P - P° 

"££ F (P-l/2)-£° F (P°-l/2)" 

= Tr (D°Q) -a [[ P -^^d X dy + ^ll PQ [ X) PQ \ v) dx dy 
JJr3 xR 3 \x-y\ 2JJ R 3 xR 3 \x-y\ 

:=^ DF (P-P°). (11) 

The energy fgrjp was introduced and studied with an exchange term by 
Chaix and Iracane in [5] (see also [B] ) . An adequate mathematical formalism 
was then provided by Bach, Barbaroux, Helffer and Siedentop [2] in the free 
case v = 0, and by Hainzl, Sere and the author in \20\ [2"T] when v ^ 0. 

The proof then reduces to finding the appropriate functional setting in 
which the minimization of the energy £g DF makes sense, providing a solution 
to 0. We quickly sketch the proof for the convenience of the reader. 



Proof. We want to give a clear mathematical meaning to the energy (jlip 
and minimize it. Let us first consider the kinetic energy term. Noticing [2 J 
that 

Q = p - p° with < P < 1 Q 2 < Q ++ - Q— 

where we have used the notation Q ee ' := P®QP®, with e, e' G {±}, we have 
(assuming all terms are well-defined), 

Tr(P>°Q) = Tt\D°\(Q ++ -Q— ) > Tr \D°\Q 2 . 

Hence the kinetic energy is nonnegative. Recalling 1 < \D°\ < y/l + A 2 on 
^Aj we also see that it is finite if and only if Q ±ziz G &i(S)a) and Q G @2(^Da)- 
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This suggests to work in the following convex subset 

K:={QtX : Q 2 <Q++ -Q—} 
where X is the Banach space 

X := {Q = Q* G 6 2 (iDA) : Q ±=t G ©i(£a)} , 
and to use the following generalized kinetic energy |20j : 
Tr P o(D°Q) :=Ti\D°\(Q++-Q—). 

Using the ultraviolet cut-off A, it was proved in \22\ Lemma 1] that the 
map Q G X i— > pq G L 2 (M 3 ) n C is continuous. Hence the energy <?b DF is 
well-defined for any state Q G KZ. 

Now, when v G C, we can complete the square and obtain the lower bound 

2 JJr3 xR 3 \x-y\ 



2 J Jr s xrz f - y| 

>-o// | , dady- (12) 

2 JJr3 xR 3 \x-y\ 

This proves that <?b DF is bounded from below. 

It is then an exercise to verify that £g DF is convex and lower semi- 
continuous on the convex set IC, hence that it possesses at least one minimizer 
Q = P — P° . It is a solution of Q with fj, = 0. Uniqueness of the density 
follows from the strict convexity of <?b DF with respect to pq. To deal with 
the case p ^ 0, one replaces D° by D° — p. □ 

The variational argument provides solutions which a priori only satisfy 
pop but one could think that they indeed have much better properties. 
As we will see in Section [H this intuition is partially wrong: solutions are 
actually quite singular. In particular Q is in general not trace-class, which 
is related to renormalization. 

The property Q ++ ,Q G &i(Sja) m (flOj) suggests to define the total 
'charge' of the system by 

q = Tr (Q ++ + Q—) :=Tr P o(Q). 

If Q is trace-class then we have Tr (Q) = Tr P o (Q) but in general Tr (Q) 
is not well-defined. Properties of the generalized trace Tr P o have been 
provided in [20]. When P = Q + P® is a projector, Trpo(Q) is always an 
integer which is indeed nothing but the relative index of the pair (P, P^ ) , 
see [201 ID- 
Varying p allows to pick the desired total charge, as we now explain. Let 
us introduce the following constrained minimization problem 

E"(q) := inf ££ DF (Q). 

Tr p o(Q)=g 

The function q i— > E l/ (q) is convex. Assume that Q = P — P® is a ground 
state for E^ DF (q). Then simple convexity arguments show that Q is also 
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Figure 2. Varying the chemical potential \i is, by convexity, 
equivalent to varying the total charge q = Tr P o (Q) of the 
system. 



a global minimizer of the free energy £g DF — //Tr po , with fi = dE v (q) / dq. 
Indeed it was shown in |17| that -Ebdf * s strictly convex on some interval 
(Qm,QM) (corresponding to ft G (—1,1)), which is also the largest interval 
on which E u (q) has ground states, see Fig. [2j Therefore varying jj, in 
(—1,1) is exactly the same as solving the minimization problem E v (q) for 
q G (<7m,<?A/)- 

Maximum ionization. The numbers q m and qM can be interpreted as the 
minimal and maximal possible ionization of the system in the presence of 
the external field v. It is important to derive bounds on these quantities, 
in order to determine for which values of the charge q the system is stable. 
The following was proved in |17| : 



Theorem 2 (Maximum Ionization [IT]). We assume that v 6 C n L 
with Z := J v > 0. 

• (Existence of neutral atoms) One has Z £ [q m ,qM]- 

• (Ionization bounds for weak external fields) For a regular ultraviolet cut- 
off, there exists constants C,C > such that the following holds: For 
a \\i/\\ c + a(l + log A) < C and Z > 0, one has 

a log A + 1/A + a\\u\\ c 
~ C 1-CalogA ^-<0< 

_ ^ 2Z + C(alogA + l/A + q|H| c ) 
- Z -' ?M - 1-CalogA • (13) 

The bound (|13p is a generalization to the BDF model of an estimate 
due to Lieb [30]. In the nonrelativistic limit (|13p reduces to Lieb's result 
Z <q M < 2Z. 

In [TTJ, the bound (fT3|) is shown by using a more regular ultraviolet cut- 
off. To be more precise, the free Dirac operator D° was replaced by an 
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operator which grows faster at infinity 

&> = #(i-±) 

and the model was settled in the whole space Sj = L 2 (M 3 ,C 4 ). 

4. Renormalization 

We have described results dealing with existence (and non existence for 
Q [lm,qM]) of solutions to Equation ([9]). All these solutions depend on 
the ultraviolet cut-off A and it is a natural question to investigate how. 
Indeed, the limit A — > oo for a fixed was shown to be very singular in |21[ 
Theorem 2] and the correct way to tackle this issue is to resort to charge 
renormalization. 

Let us start by recalling the spirit of renormalization. A physical theory 
usually aims at predicting physical observables in terms of the parameters 
in the model. Sometimes, interesting quantities are divergent and it is nec- 
essary to introduce cut-offs. In our case the parameters are the coupling 
constant a = e 2 , the cut-off A, the chemical potential fj, and the external 
density v. For simplicity we will take \x = and we will not emphasize 
the dependence in v in our notation. In our system of units the mass of 
the electron is m = 1. Predicted physical quantities are functions F(a,A). 
The charge e (or equivalently its square, the coupling constant a) is also a 
physical observable and renormalization occurs when the value predicted by 
the theory is different from its 'bare' value: 

a p h = a ph (a, A) / a. (14) 

In this case the parameter a is not observable in contrast with a p h = 
a p h(m, a, A) which has to be set equal to its experimental value. The rela- 
tion ()14p has to be inverted, in order to express the bare parameter in terms 
of the physical one: 

a = a(a ph ,A). (15) 

This allows to express any observable quantity F as a function F of the 
physical parameters and the cut-off A: 

F(a ph ,A) = F(a(a ph ,A), A) (16) 

A possible definition of renormalizability is that all such observable quanti- 
ties have a limit when A — > oo, for fixed a p h- 

Important difficulties can be encountered when trying to complete this 
program. For instance the physical quantity a p h might be a nonexplicit 
function of a. The corresponding formulas can then only be inverted per- 
turbatively to any order, as is the case in QED [i~2"l H], [27] . 

Nonperturbative charge renormalization formula. For the model pre- 
sented in this article it is fortunate that there is an explicit and nonpertur- 
bative relation between a p h and a, as expressed in the following result: 

Theorem 3 (Nonperturbative charge renormalization formula [17]). As- 
sume that a > ; A > and /i £ (—1,1) are given and let P be a solution 
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of © as given by Theorem [0 If v £ C n L 1 



, i/ien p P _po G L 1 



Pp-pO 



[ u-Tr P o{P-P° 
1 + <xB A 



where 



B A 



7T 



z 4 /3 



1 



3vr 



>o g A-A + ^ + 0( V A>) 



(17) 



(18) 



Note that, except in the neutral case Tr P o(P — P_) = Z, (fTTj) implies 

that the solution Q = P — P® found in Theorem [1] cannot be trace-class. 
If Q were trace-class, we would have Tr P o(Q) = Tx (Q) = f R3 pq which 
contradicts (fTT|) . 

The previous result is interpreted as follows. Assume that we put a nu- 
cleus of charge Z in the vacuum, and let P be the corresponding Dirac's 
polarized vacuum (that is we take p = in ©). When a |H| C is small 
enough n it was proved in [21] that it holds \\P - P°|| < 1. This itself 
implies that the relative index vanishes, Tr P o (P — P°_) = 0, hence 



Pp-pO 



Z 



1 + aB A 

In reality we never measure the charge of the nucleus alone, but we always 
also observe the corresponding vacuum polarization. Hence the physical 
coupling constant is given by the renormalization formula 



Ctph 



a 



1 + aB A 



a 



1 - a P h#A 



(19) 



In our theory we must fix a p ^ and not a. Using f)19|) we can express any 
physical quantity in terms of a p h and A only. 

Unfortunately it holds a p ^B\ < 1 hence it makes no sense to take A — > oo 
while keeping a p h fixed (this is the so-called Landau pole [28]) and one has 
to look for a weaker definition of renormalizability. The cut-off A which was 
first introduced as a mathematical trick to regularize the model has actually 
a physical meaning. A natural scale occurs beyond which the model does not 
make sense. Fortunately, this corresponds to momenta of the order e 37r / 2 "p h , 
a huge number for a p h ~ 1/137. 

Asymptotic renormalization. It is convenient to define a renormalized 
density p p h by the relation [21] 

a p hP P h = a{v - p P _po) (20) 

in such a way that D = D° — a p h/> p h * \%\~ ■ This procedure is similar to 
wavefunction renormalization. By uniqueness of Pp^po we can see p p h as 
a function of a p h, v, p and A (or n). For the sake of clarity we do not 
emphasize the dependence in v and we take p = (this means that we 



^By scaling we can keep f R3 v fixed and choose a ||f[ c as small as we want. 
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consider the vacuum polarization in the presence of the nucleus, without 
any real electron). The self-consistent equation for /? p h was derived in |21j . 
It is explained in [TSj that one can expand p p h = /9 p h(a p h>A) as follows: 

oo 

Pph(Qph)A) = y^(o;ph) n ^n,A (21) 

n=0 

where {v n ,\}n C L 2 (M 3 ) n C is a sequence depending only on the external 
density v and the cut-off A. This sequence is defined by an explicit induction 
formula which is detailed in [18J and that we do not write here for short- 
ness. The series (|21 j) has a positive radius of convergence, which is however 
believed to shrink to zero when A — > oo. 

Assuming v decays fast enough (see condition ([22j) ). it is proved in [18] 
that for any fixed n, the limit v n ,A —> v n exists in L 2 (IR 3 ) f)C. This is what is 
usually meant by renormalizability in QED: each term of the perturbation 
series in powers of the physical a p h has a limit when the cut-off is removed. 
The sequence {u n } n is the one which is calculated in practice [4"1 [19 ], [TE [ [T4"] . 
One has for instance Uq = v and 

^ * M-l = _L / dtit 1 -!) 1 ' 2 
3vr J x 

the Uehling potential I35| . All the others v n can be calculated by in- 
duction in terms of vq, v n -\. An explicit recursion relation is provided in 

The next natural question is to understand the link between the well- 
defined, cut-off dependent, series (f2Tj) and the formal series S^Lo( a ph) ni/ V 
Recall that a p h-BA < 1 by construction, so it is in principle not allowed to 
take the limit A — > oo while keeping a p h fixed. 

It is more convenient to change variables and take as new parameters 
a p h and k = a p ^B\, with the additional constraint that < k < 1. The 
new parameter k is now independent of a p h and we ask ourselves whether 
predicted physical quantities will depend very much on the chosen value of 
< k < 1. The purpose of [18] was precisely to prove that the asymptotics 
of any physical quantity in the regime a p h <C 1 is actually independent of k 
to any order in a p h, which is what was called asymptotic renormalizability. 
Note that fixing k £ (0, 1) amounts to take an a p h-dependent cut-off A ~ 

(j e 3nK/2a ph ^ 

Theorem 4 (Asymptotic renormalization of the nuclear density [IB])- Con- 
sider a function v G L 2 (M 3 ) n C such that 

[ log(l + |fc|) 2JV+2 |P(fc)| 2 dfc < oo (22) 

for some integer N. Let p p h(a P h> K ) be the unique physical density defined 
by (|20p with \i = 0, a p h > and < k < 1, corresponding to the bare 
coupling constant a = (1 — K) _1 a p h and the ultraviolet cut- off A such that 
Ba = K/a p h- 



t 2 + t 4 



\x - y\ 
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Then, for every < e < 1, there exist two constants C(N,e,v) and 
a{N,e,u), depending only on N, e and v, such that one has 



for all < a p h < a{N, e, u) and all e < k < 1 — e. 

The interpretation of Theorem[3]is that the renormalized density p p h(ctph, K ) 
is asymptotically given by the formal series Y2n>o( a ph) n v n , uniformly in the 
renormalization parameter k in the range e < k < 1 — e. For a very large 
range of cut-offs, 



the result is independent of A for small a p h- This formulation of renormal- 
izability is more precise than the requirement that each u n \ converges. It 
also leads to the formal perturbation series in a very natural way. 

It was argued by Dyson in [T3] that the perturbation series X^n>o( a ph) n ^ 
it is probably divergent, but there is no mathematical proof so far. In [18j . 
some properties of the sequence {f n } were derived. 



We have presented a mean- field theory for electrons in atoms and molecules, 
which describes at the same time the self-consistent behavior of Dirac's po- 
larized vacuum. The model can be deduced from Quantum Electrodynamics 
by restricting to Hartree-Fock states and neglecting photons in the Coulomb 
gauge. 

Existence of ground states could be established, with or without a charge 
constraint. The so-obtained states are rather singular, in particular they 
yield a perturbation Q of the free vacuum P® which is in general not trace- 
class but still has pq G L 1 (1R 3 ). This technical issue is at the origin of charge 
renormalization. 

The formula linking the physical coupling constant a p h and the bare a 
is explicit and exhibits a Landau pole, rendering impossible to remove the 
ultraviolet cut-off A while keeping a p h fixed. Nevertheless in a regime where 
a ph 1 and A S> 1 such that k = (2/37r)a p h log A stays bounded, the 
asymptotics is found to be independent of the value of k, to any order in 
the physical coupling constant a p h- The terms of the asymptotic expansion 
are the ones which are computed in practice. The first order term induces 
the famous Uehling potential. 

The model which we have presented in this paper is probably not quan- 
titative but it already possesses several of the qualitative properties of full 
Quantum Electrodynamics, with the advantage that they can be studied in 
a fully rigorous manner. A more quantitative model would include photons, 
for instance via an additional self-consistent classical magnetic field, as is 
done in Relativistic Density Functional Theory. 




(23) 



£< ie 3e7r/2a ph < ^ < (J 2e 3{l-e)n /2a ph 



Conclusion 
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